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■ Abstract. Over the last few years, high-quality X-ray imaging and spectroscopic 

data from Chandra and XMM-Newton have added greatly to the understanding of 
the physics of radio jets. Here we describe the current state of knowledge with an 
emphasis on the underlying physics used to interpret multiwavelength data in terms 
of physical parameters. 
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^ ■ 1 Introduction 

(N 

o 

Jets in active galaxies emit over a wide range of energies from the radio to the 
7-ray. Synchrotron radiation and inverse-Compton scattering are the two main 
■ radiation processes, with their relative importance depending on observing 

^ <— | frequency and location within the jet. The thermally-emitting medium into 

Qh' which the jets propagate plays a major role in the properties of the flow and 

the appearance of the jets. 
5^ ■ In this chapter, we concentrate primarily on the interpretation of obser- 

C/j , vations of resolved jet emission. Much new information about jets has been 

gained over the past few years. The high spatial resolution of Chandra has 
been key to the study of X-ray jets, and the large throughput of XMM-Newton 
has assisted studies of the X-ray-emitting environments. While historically the 
search for optical jet emission has been carried out using ground-based tele- 
scopes, a major problem has often been one of low contrast with light from 
the host galaxy. The sharp focus of the Hubble Space Telescope (HST) helps 
to overcome this difficulty, and the HST is playing a particularly important 
role in optical polarization studies. The recent opening of the Very Large 
Array (VLA) data archive has greatly assisted multiwavelength studies of ra- 
dio jets. It, together with the Australian Telescope Compact Array (ATCA) 
in the southern hemisphere, allows multifrequency radio mapping at angular 
resolutions well matched to Chandra and XMM-Newton. The combination of 
data from the radio to the X-ray has been a key element in advancing our 
understanding of jets. 
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Active-galaxy jets provide an exciting laboratory where the interplay of 
relativistic effects, plasma physics, and radiation mechanisms can be seen. Our 
purpose in this chapter is to provide an introduction to the field and some 
pointers to recent research. 

2 Radiative processes 

The non-thermal mechanisms of synchrotron radiation and inverse Compton 
scattering are described extensively in several places, e.g., |14H49ll5Tll96lllllj . 
and descriptions of thermal radiation can be found, e.g., in |81llllll|ll4U123] . 
In the following sections we give equations important in the interpretation of 
jet emission in useful forms that are normally independent of the system of 
units. 

2.1 Synchrotron radiation 

The rate of loss of energy of an electron (or positron) of energy E is given by 

- dE/dt = 2a T cj 2 p 2 u B sin 2 a (1) 

where <tt is the Thomson cross section, c is the speed of light, 7 is the Lorentz 
factor of the electron (= E/m e c 2 ), (3 is the speed of the electron in units of c, 
Ub is the energy density in the magnetic field, and a is the pitch angle between 
the direction of motion and the magnetic field. Averaging over isotropic pitch 
angles, P(a)da — (1/2) sin a da, 

- dE/dt = (4/3)cr T c7 2 /3 2 w B (2) 
The radiative lifetime of the electrons is usually calculated as 

Trad = E /{-dE/dt) (3) 

which, for energy losses proportional to energy squared, as here, is the time for 
any given electron to lose half of its energy. High-energy electrons, responsible 
for high-energy radiation, lose their energy fastest. 

The spectral distribution function of synchrotron radiation emitted by 
monoenergetic electrons of Lorentz factor 7 is rather broad. It is usual to 
define the critical frequency, 

v c = (3/2) 7 2 ^ g sina (4) 

where v g is the non-relativistic electron gyrofrequency, which is proportional 
to the magnetic field strength, B. Written in SI units, v g = eB/2?:m c ss 
30£> GHz, where B is in units of Tesla. As a rough approximation, something 
close to v c can be used as the frequency of emission, but for the full distribution 



Multiwavelength evidence of the physical processes in radio jets 



3 




Fig. 1. The spectral distribution function of synchrotron radiation from a monoen- 
ergetic electron of Lorentz factor 7, where u c is defined in JIJ. The spectrum is 
plotted linearly on the left and logarithmically on the right 



function it is convenient to define X = v jv Cl and then the spectral distribution 
function depends on frequency through 

/>oo 

F(u,u c ) = X K 5/3 (()d( (5) 
Jx 

where K 5 / 3 is the modified Bessel function of order 5/3. The spectrum peaks 
at 0.29f c , as shown in Fig. ^ The synchrotron spectrum from a distribu- 
tion of electrons with some maximum Lorentz factor, 7 ma x) will therefore fall 
exponentially at high frequencies. The gradual turn-down at low frequencies, 
oc v 1 / 3 , is not expected to be seen in practice. Instead, for a homogeneous emit- 
ting region the low-energy fall off will have a slope oc v 2 b due to the source 
having become optically thick and suffering from synchrotron self-absorption. 

The luminosity per unit frequency for a number spectrum of electrons 
N{"i)d~f is 

c m c r c z^ g sin a / F(v,v c ) Nij)^ (6) 

where r e is the classical electron radius, and this form is valid in all systems of 
units. For a power-law number distribution of electrons N("f) dj = kj~ p dj, an 
analytical result can be found for frequencies satisfying 7mi n ^g ^ v ^ Tmax^g 



(7) 

where r is the gamma function, and, again this form is valid in all systems 
of units. If the pitch angle distribution is isotropic, the term in sin a may be 
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Fig. 2. The spectral distribution function of inverse Compton radiation from a mo- 
noenergetic electron of Lorentz factor 7 in an isotropic radiation field of frequency 
v a in the Thomson limit JHJ. The spectrum is plotted linearly on the left and loga- 
rithmically on the right 



replaced by 0rf[(p + 5)/4]/2f[(p + 7)/4]. The result for more complicated 
electron spectra must be computed numerically using 

It is usual to use the symbol a for the negative exponent of the power-law 
radiation spectrum [L v oc v~ a ), where a = (p— l)/2. Particle acceleration by 
ultra-relativistic shocks, likely to be important in jets, produces p w 2.2 — 2.3 
1 , and p will be steeper in regions where energy losses are important. 

2.2 Inverse Compton scattering 

Relativistic electrons lose energy by scattering photons to higher energy in 
a process analogous to synchrotron radiation, where virtual photons are re- 
placed with real ones. Providing an electron of mass m e and Lorentz factor 7 
is scattering a photon of low enough energy such that 

"fhv Q <C m c c 2 (8) 

the Thomson cross section, ot is applicable, and the rate of loss of energy of 
the electron in an isotropic radiation field of total energy density u ra( j is 

- dE/dt = (4/3)a T C7 2 /3 2 Urad (9) 

where other symbols are as defined in Sect. 12. ll The radiative lifetime can be 
calculated using ©. Where electrons are losing energy both by synchrotron 
and inverse Compton emission, the ratio of total luminosity in the two emis- 
sions is simply L iC /L syn = u ra d/uB- 

An exact result ^31 exists for the spectral distribution function for a 
monoenergetic electron of Lorentz factor 7 scattering an isotropic radiation 
field of photons of frequency v a up to frequency v. It is convenient to define 
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Fig. 3. The combinations of electron Lorentz factors and photon frequencies which 
will produce photons of typical optical, X-ray and 7-ray energies via inverse Compton 
scattering of an isotropic photon distribution in the Thomson limit is shown on the 
left. The Thomson limit only holds if fhuo/m e c 2 <C 1. The right-hand plot shows 
this is satisfied for optical and X-ray photons, but for 7-ray and higher energies the 
Klein-Nishina cross-section begins to become important 



X = v/A^Vo, where (I/47 2 ) < X < 1, and then the spectral distribution 
function is proportional to 

F(v,u ,j) =Xf(X) = X(l + X -2X 2 + 2X\nX) (10) 

This is plotted in Fig. [21 The luminosity per unit frequency for a number spec- 
trum of electrons N(^)drf and a spectral number per unit volume of photons 
n(y )dv is 

L„ = 3ft,ecr T J J n(v ) N(i) F(v, v , 7) dv drf (11) 

The relationship between the electron power-law index, p, and the spectral in- 
dex of inverse Compton radiation, a, is the same as for synchrotron radiation. 
Since the mean value of X, 

Xf(X)dX/ [ f{X)dX = 1/3 (12) 



the mean frequency of photons scattered in an isotropic radiation field is given 

by 

v = (4/3)7 2 ^ (13) 

Fig. |3 plots combinations of 7 and v a that will produce optical, X-ray, 
and 7-ray photons according to (|13fl . and shows that in most situations (JSJ) 
is applicable. As © becomes violated, the Klein-Nishina cross section must 
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be used in place of <7t, as in these situations the electron begins to lose a 
significant fraction of its energy in a single scattering, e.g., [T4"| . 

The photon fields available to the electrons in jets are the synchrotron 
radiation they have produced, in which case the emission is known as syn- 
chrotron self-Compton (SSC) radiation, the cosmic microwave background 
(CMB) which scales with redshift as (1 + z) A , and photons from the active- 
galaxy nucleus. In practice the photon fields are rarely expected to be isotropic 
when the effects of the geometry and bulk speed of the jet plasma are taken 
into account. 



2.3 Thermal radiation and energy loss 

For hot X-ray emitting gas containing ions of charge Zj, and where the cooling 
in line radiation is unimportant, the X-ray emissivity, £ „ depends on temper- 
ature, T, and electron and ion number densities, 

327r / 2tt\ 1/2 n 23 ( m c c 2 \ 1/2 _ hv/kT M iS 

V = ~\J) Z C9snin c m c c 2 r 3 c [-j^rj * 1 (14) 

for all systems of units, where r e is the classical electron radius. gg(y,T, Z), 
the free-free Gaunt factor which accounts for quantum-mechanical effects, is 
of order unity and a weak function of frequency. An approximate form |1 1 1| 
in the case of most X-ray interest, (kT/eV) > 13.6Z 2 and Z\ < 2, is 

hv>kT 

<m = \ ,v;::vum hv<kT ( l5 ) 

where the constant £ = 1.781. For heavier ions more complicated forms must 
be used, and since the heavier ions contribute disproportionately to the X-ray 
output (because of the Zf factor), the calculation of the correct Gaunt factor 
becomes a computational issue. 
In SI units, (|14|) becomes 




£ v = 6.8xlQ- 51 Z?T- 1/2 (n c /m- 3 )(n i /m- 3 )g s e- hu/kT W nT 3 Hz" 1 (16) 

At temperatures below ~ 2 keV, line radiation cannot be ignored, and 
the plasma models incorporated into X-ray spectral-fitting programs such as 
XSPEC 0] can be used to find £ v . Examples for plasma of kT = 0.5 keV and 
5.0 keV are shown in Fig. 0] It is common to define the rate of loss of energy 
per unit volume per unit frequency per n p n e as A(v), where n p is the proton 
number density. In XSPEC, A(v) is normalized by (1 + z) 2 J n p n c dV '/ 'AnD 2 ^ in 
fitting to data, where V is volume and Z?l is luminosity distance. This allows 
n p to be determined assuming some geometry for the source, as illustrated in 
Sect.il 

The total energy-loss rate per unit mass is given by 
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Fig. 4. Emissivity distributions of gas of kT — 0.5 keV (left) and kT = 5 keV (right) 
with normal cosmic abundances from the APEC model in XSPEC. The line emission, 
superimposed on the thermal bremsstrahlung continuum, is a more prominent chan- 
nel for cooling in cool gas. The resolution of these model spectra is higher than that 
of current X-ray detectors 



e= xj^ (17) 

m H n P 

where X is the mass fraction in hydrogen, which is 0.74 for normal cosmic 
(i.e., solar) abundances. 

If the energy loss is sufficiently rapid, the radiating gas will fall towards 
the centre of the gravitational potential well as its pressure support fails. The 
rate of infall is therefore determined by whether the characteristic timescale 
for radiative energy loss is greater than or less than the timescale for energy 
input, perhaps from the dissipation of energy by the radio jet. 

If the timescale for energy loss is defined by the logarithmic rate of change 
of temperature in the comoving fluid frame, 

where D/Dt is the convective derivative (d/dt + v ■ V), which describes the 
rate of change in a frame moving at velocity v with the gas, then the value of 
t coo i can be obtained from the entropy equation 
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kT (keV) 

Fig. 5. Curves of r coo in p (with n p in units of 10 3 m -3 or 10 -3 cm -3 ) as a function 
of kT for three metal abundances with respect to solar: 0.1, 0.5, 1.0 



where s(T, P) is the entropy per unit mass of the gas and P is the pressure. 
For isobaric (P = constant) losses in a monatomic gas, 

5 kT 

Tcool = 77 20 

Curves of r coo in p as a function of kT, calculated using the APEC plasma 
model in XSPEC, are shown in Fig. At high temperatures, r coo i cx T 1 / 2 , but 
at T < 2 keV line emission from common metal species in the gas becomes im- 
portant, and so e and r coo i become strong functions of kT. The offset between 
the curves at large kT arises from the increased density of highly-charged ions 
in the radiating plasma as the metal abundance increases. 



3 Effects of bulk relativistic motion 

The first evidence for bulk relativistic motion in jets was on pc scales. In 
the early 1970s, the technique of Very Long Baseline Interferometry (VLBI) 
became established, and the first detections of apparent superluminal mo- 
tion of small-scale radio components were being made (see [57] for a short 
review). After some years of debate, it became widely accepted that the true 
explanation of the effect was as proposed earlier, |1U7| . and that the appar- 
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ently superluminal velocities resulted from components moving towards the 
observer at bulk relativistic speeds. 

From a consideration of light travel time, it is straightforward to show 
that a source moving with a relativistic speed (3c, bulk Lorentz factor r (= 
1 / y 1 — (3 2 ) , at angle 6 towards the observer, has an apparent transverse ve- 
locity of 



Uapp =/3csin0/(l -ftcasd) (21) 

which takes on a maximum value of F(3c when sin# = 1/-T. Hence large (3 
(> l/v2) and small 6 can easily result in w ap p > c. 
The bulk relativistic Doppler factor is defined as 



5 = l/T(l-/3cos(9) (22) 

and takes on a large value for large (3, small 8. For a spherical blob of emission, 
the apparent spectral luminosity, L v , is increased by <5 3+Q , where a is the spec- 
tral index fSect. |2~T|) . although in most steady jets 5 2+a is more appropriate 

unsung. 

A variable source will show intensity changes on timescales shorter by S 
than the true value, and beaming effects can have a dramatic effect on model 
parameters. For example, it was shown that the predicted X-ray flux density 
from SSC in a compact, spherical, variable radio source must be reduced 
by a factor of (5 2 ( 3q + 5 ) if variability is used to measure the size, and the 
observed self-absorption turn-over frequency is used to estimate the magnetic- 
field strength [SSI EH] • 

Relativistic beaming has a major effect on a jet's appearance, and by the 
1980s the first models to unify classes of active galaxies through jet orientation 
were proposed, e.g., |94M116| . and subsequently further developed to account 
for multiwavelength properties, e.g., |rJl 11151 fl30| . It is common to assume 
that a radio source has two oppositely-directed jets that are intrinsically the 
same, although there is doubt that this need be the case either from jet- 
production considerations or asymmetries in the environment through which 
the jets propagate. Under the assumption of intrinsic similarity, the brightness 
ratio between the approaching and receding jets is given by 



1 + /3cos6>, 

This constraint on (3 and is often used in conjunction with a constraint 
on u app to estimate separately the jet speed and the angle to the line of sight, 
e.g., |fl8lll47j . The core dominance, R c< \, defined as the flux density of beamed 
emission divided by that of unbeamed extended emission, is also commonly 
used as an indicator of the angle of a source to the line of sight. When applied 
to a twin-jet source, e.g., 
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-Red, 



(cob 8 = 0) 



(1 - /5cos6»)-( Q + 2 ) + (1 + f3 cos 6)- ("-V 



(24) 



2 



Bulk rclativistic motion towards the observer has an important effect on 
inverse-Compton scattering in the case of an external photon field (i.e., not 
SSC) 34 . In its rest frame an electron now sees a directionally enhanced 
radiation field and the rate of scatterings increases. If the photon field is 
isotropic in the observer's frame (such as the CMB), the increase in scattering 
rate is by a factor of 5. Since lower-energy electrons (of which there are more) 
are needed to scatter photons of a given rest-frame frequency to the observed 
frequency, there is a further enhancement by 5 a . This means that the observed 
luminosity at a given frequency from inverse-Compton scattering of the CMB 
will be increased relative to the synchrotron luminosity by 5 1+a , assuming 
that the magnetic field and electron spectrum are held constant. Although 
this in itself is not a large factor compared with ~ 6 2+a (see above) that 
arises from any radiation mechanism and applies both to the synchrotron 
and inverse-Compton emission, it can be sufficient to raise the ratio of X-ray 
to radio emission enough to fit observations of quasar jets (Sect. 16. 3[) . For a 
source radiating at minimum energy fSect. RTT|l . while the X-ray emission from 
inverse-Compton scattering of CMB photons is increased, the X-ray emission 
from SSC is reduced. To increase this ratio via the SSC mechanism requires 
taking the source out of minimum energy through adding relativistic electrons. 
A beamed source has intrinsically lower total energy content, and, despite the 
higher jet speed increasing the kinetic energy, it may be in a condition of 
minimum total (radiative and kinetic) power |47) . 

4 The external medium 

Often the bulk of the X-ray emission from an active galaxy is not jet re- 
lated, but rather arises from the hot interstellar, intergalactic, or intracluster 
medium around the jet. If we assume that this gas is almost at rest in the 
local gravitational potential well, with little transfer of energy or matter to or 
from it, then it will be in a state close to hydrostatic equilibrium. 

Under these circumstances, the gas will take up a density and temperature 
profile which depend on the distribution of mass and the thermal history of 
the gas. The principal governing equation is that of hydrostatic equilibrium, 
with 



where is the gravitational potential and p is the total mass density of the 
gas. We may assume that the gas has the ideal gas equation of state 



VP 



(25) 



P = 



pkT n p kT 



(26) 
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to relate P and p, where p is the mass per particle in units of mu, given by 



3.. 1 x 1 



l-i = \^2X + -Y + -Z\ = 0.60 (27) 

for a gas with solar abundances (X = 0.74, Y = 0.25, Z = 0.01). 
If the mass distribution is spherically symmetrical, then 



d£ = GMtotjr) 
dr r 2 



(28) 



where M to t(r) is the total mass within radius r, and so the mass density and 
temperature profiles of the gas are related to the distribution of total mass by 

kTr f d In p d In T 



pm p 



GM tot (r)=-^(^ + — ] . (29) 



A consistent solution of this equation is obtained by assuming that the gas 
is isothermal and that the mass distribution follows 



Mtot (r ) = 2 M c (30) 
r c (r A + r~) 

where r c , the core radius, defines the characteristic scale of the mass distri- 
bution and M c is the mass within r c . The density of the gas then follows the 
so-called isothermal (3 model 

P = Po[l + ^J (31) 

where (3 is a constant which determines the shape of the gas distribution, and 
depends on the ratio of a characteristic gravitational potential to the thermal 
energy per unit mass in the gas 

d kl r c 

The alternative derivation of (|31|l given by brings out the interpreta- 
tion of P in terms of the relative scale heights or, equivalently, the ratio of 
energy per unit mass in gas and dark matter. p , r c , and [3 in (|31fl are of- 
ten fitted to provide convenient measures of the gas distribution's mass, scale 
and shape without considering the detailed properties of the underlying mass 
distribution (|30fl . 

It should be noted that the physical consistency of this much-used model 
for the gas distribution depends on radial symmetry (a simple distortion into 
an ellipsoidal model for the gas density p would imply a mass distribution 
which is not necessarily positive everywhere), and on gas at different heights 
in the atmosphere having come to the same temperature without having neces- 
sarily followed the same thermal history. This is problematic given the various 
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origins of gas — some will have fallen into the cluster, some will have been 
ejected from stars, and some may have been moved by the jet. The solu- 
tion 1|31|) is not unique in the sense that gas with a different thermal history 
might follow a significantly different density distribution. Thus, for example, 
if the gas has the same specific entropy at all heights, and so is marginally 
convectively unstable, with 



P oc p 1 (33) 

with 7 being the usual specific heat ratio (5/3 for a monatomic gas), the mass 
distribution (|30|l would lead to a gas density distribution 



P'-'^pT U-/3Aln l + (34) 



where /3a is a structure constant with a similar meaning to /3, and the relation 
only applies within the radius at which p — ► 0. 

A similar procedure for an NFW mass profile and an isothermal gas 
leads to a density distribution 

/ \ ar s /r 

P = p s 2- a fl + f J (35) 
where the new structure constant is 

1 pmn GM S 

a= hT2^ri^^r (36) 

(r s is the scale of the NFW model and M s is the mass within radius r s ). It 
can be seen that in this solution p — > oo as r — > 0, since the NFW density 
profile has a cusp at r = 0. 

The total masses for either the NFW profile or (|30|l diverge as r — > oo, 
so both profiles must be truncated at some radius, such as r2ocb the radius 
at which the mean enclosed mass density is 200 x the critical density of the 
Universe, p cr it at redshift z. Care should be taken in using the solutions for 
p to ensure that the fraction of mass in gas does not exceed the cosmological 
bound (~ 17%). 

The run of density and temperature in the atmosphere around a jet are 
measured from the X-ray image and spectrum, where a density model for the 
gas of the form of l|31(l or l|35[) is fitted to the X-ray surface brightness. The 
X-ray surface brightness at frequency v at a point offset by r in projected 
distance from the centre of a spherical gas distribution is (in energy per unit 
time per unit solid angle per unit frequency) 

Z„{r) d n = { l + z) t n ^p dV (37) 
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where A v is the emissivity of the gas at frequency is, and is the luminosity 
distance of the gas. The (1 + z) factor takes account of the redshifting of time 
and frequency in the definition of S v . The volume element in the integral is 

dV = dlD\dQ (38) 

where dl is an element of distance along the line of sight, Da is the angular 
diameter distance, and dfi is the element of solid angle, so that the surface 
brightness becomes 

, „ f n r n n A,,{T)dl 

s ^= j 4 ; ( i + v • ^ 

For an isothermal gas with a /3-model density distribution 

„2 , 7 2 \ -3/5 



4n(l + z) 3 J V r c ) 



dl (40) 



where the integral is taken along the line of sight with limits of ±oo or to 
some cut-off radius. In the former case the integral is simple, and the X-ray 
surface brightness can be written JT] 



Sv{r) = 47 r(l + z)3 + TW) < 41 > 

It is also useful to write down the observed flux density, S u , of a gas distribu- 
tion out to angle 9 from the centre, since this is the quantity generally fitted 
in X-ray spectral analyses. This is 



o 



S u (9) = I 2n0de£ v (r) (42) 

02 X i-3/J\ 

1 + 7^ (43) 



A V (T) n c0 n p0 8 3 C D A ^ r{i(3 - f) / 
(1 + z) 3 4 r(3/3) ^ 

where C is the angular equivalent of the core radius r c . Now, XSPEC (Sect. l2~3*|) 
calculates a normalization, A/", for any thermal gas model, where 

Jn c n p dV _ (l + z)S v 
~ 4^1(1 + zf ~ A v [ > 

and this is related to the parameters of the isothermal /3 model by 

n c0 n p0 e 3 c D L VI r(3f3-l) ( _ f *\ ^ 3/3N 



M = cu pu c ^ s± l _ l + _ (45) 

a + z)4 4 nm [ { + 61) ) 1 ° j 

The fitting process in XSPEC takes account of the form of A v and the 
conversion from energy units (£„) to the count rate in energy bins used by 
X-ray detectors. 
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This expression shows how to use the normalization Af found by XSPEC, 
or some other code, from the X-ray spectrum within angular radius 9 of the 
centre, and (3 and r c as fitted from the radial profile of the X-ray emission, to 
measure the central electron and proton densities which are related by 

± = ^ ~ 1.18 (46) 

The proton density found in this way, and the gas temperature measured 
from the spectrum, provide crucial information on the gas environments of 
radio jets and hence the external forces on the jet. The separation of X- 
rays from the non-thermal processes associated with the radio source and the 
thermal processes associated with the environment is therefore an important 
part of establishing the physics of a jet. This is particularly true for low power 
(FRI) jets, which are in direct contact with the external medium fSec. \b.'2\ . 
but is also important for establishing the balance of static and ram-pressures 
in powerful sources. 



5 Simple radio-jet models 

The FRI and FRII classifications, into which extragalactic radio sources were 
separated by remain in wide usage today. Although based on radio mor- 
phology, 02 found a relatively sharp division in radio luminosity, with most 
sources of total 178-MHz luminosity below 2 x 10 25 W Hz -1 sr _1 being of FRI 
classification, and the others of FRII. Some dependence of the location of the 
FRI/FRII boundary on the optical magnitude of the host galaxy has also been 
found [HH] , suggesting the effects of environment or the mass of the originating 
black hole are important in determining the type of radio galaxy formed by 
an active nucleus. However, the primary cause of the structural differences is 
believed to be due to the speed of the primary jet fluid, with the beams pro- 
ducing FRII sources being supersonic with respect to the ambient medium, 
whereas FRI jets are seen as the result of turbulent transonic or subsonic 
flows. The following sub-sections describe the observational consequences. 

5.1 High-power FRII jets 

The standard model for the expansion of a powerful radio source powered 
by a jet which is supersonic with respect to the X-ray-emitting interstellar 
medium (ISM) is illustrated in Fig. [fj] The jet terminates at the beam head 
(in a feature identified as the radio hotspot) where the jet fluid passes through 
a strong shock to inflate a cocoon of radio-emitting plasma. The energy and 
momentum fluxes in the flow are normally expected to be sufficient to drive 
a bow shock into the ambient medium ahead of the jet termination shock. In 
the rest frame of the bow shock, ambient gas is heated as it crosses the shock 
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Fig. 6. In the standard model for powerful radio sources, a supersonic jet (dark 
grey) terminates at the beam head, producing a radio hotspot. Provided the shocked 
radio-emitting fluid forming the radio lobe (light grey) has enough internal energy 
or momentum density to drive a leading bow shock, ambient X-ray-emitting gas will 
be heated as it crosses the bow shock to fill the medium-grey region 




Fig. 7. 8.4 GHz VLA A and B array radio image of an example high-power FRII 
source, the radio galaxy 3C 220.1. The resolution is 0.3 x0.2 arcsec, and the image is 
from |145| . The angular length of the source is about 28 arcsec which, at the galaxy's 
redshift of z = 0.61, corresponds to 187 kpc (Ho = 70 km s _1 Mpc -1 , Q m — 0.3, 
Q A = 0.7) 
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to fill a region surrounding the lobe of radio-emitting plasma. Observationally 
on kpc scales we see the radio emission from well-collimated jets feeding edge- 
brightened lobes. An example is shown in Fig. [7| 
The sound speed in gas of temperature T is 



cs = Jl™ (47) 

where 7 is the ratio of specific heats (7 = 5/3), k is the Boltzmann constant, 
mn is the mass of the hydrogen atom, and /imn = 0.6mn is the mass per parti- 
cle. Under these conditions, c s w 516(fcT/keV) 1 / 2 km s" 1 w 0.54(fcT/keV) 1 / 2 
kpc Myr -1 . 

The Mach number of the speed of advance, u a dvj of the bow shock into the 
ambient medium is M. = v a dv/c s , which in convenient units can be expressed 
as 

M ra 580(w ad v/c)(/cT/keV)- 1/2 (48) 

where c is the speed of light. The advance of the bow shock is likely to be slow 
with respect to the bulk speed of the jet. However, the jets need not have bulk 
relativistic motion (Sect. for the general description of an FRII source to 
hold. 

The state of the X-ray gas close to a radio galaxy with a leading bow shock 
should reflect the source dynamics, and in addition to the X-ray emission of 
the ambient medium (Sect.0J), we may also expect to see hotter shocked gas 
surrounding the radio lobe (Fig. EJ) - In a simple application of the Rankinc- 
Hugoniot conditions for a strong shock |123| . the pressure, density, and tem- 
perature, respectively, in the unshocked (subscript 1) and shocked (subscript 
2) regions at the head of the bow shock are related by 

Pi/ Pi = (5M 2 - l)/4 (49) 
P2/P1 = AM 2 /{M 2 + 3) (50) 

T 2 /Tx = (5M 2 - 1){M 2 + 3)/16M 2 (51) 

for a monatomic gas. 

Using the above equations, and converting density into an emissivity using 
XSPEC, as discussed in Sect. 12.31 we find that in the energy band 0.8-2 keV, 
where Chandra and XMM-Newton are most sensitive, for an Ai — 4 shock, 
the X-ray emissivity contrast between shocked and unshocked gas is a factor 
of 3 higher if the ambient gas is at a galaxy temperature of ~ 0.3 keV than if 
the external medium has a cluster temperature of ~ 4 keV. 

Complications apply in reality, and in practice these are difficult to treat 
even with data from observatories as powerful as Chandra and XMM-Newton. 
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Fig. 8. Radio images of an example low-power FRI source, the radio galaxy 3C 31. 
The left image shows the large-scale structure. The map |8U| has a resolution of 
29 x 52 arcsec and is from a 608 MHz observations using the Westerbork Synthesis 
Radio Telescope. The angular length is about 40 arcmin which, at the galaxy's 
redshift of z = 0.0169, corresponds to about 830 kpc (H = 70 km s _1 Mpc -1 ). The 
right image shows an 8.4-GHz VLA map of the the inner ~ 18 kpc where the jets 
are relatively straight. The resolution is 0.25 arcsec, and the data are from |78| 

Firstly, there is observational evidence that in supernova remnants the post- 
shock electrons are cooler than the ions [23 111)6] . Secondly, the simple 
Rankine-Hugoniot equations do not take into account the fact that the bow 
shock around a lobe is oblique away from its head, with a consequent change 
in the jump conditions and the emissivity contrast |135| . However, in Sect. I6.ijl 
we describe how the above equations can also be applied to an overpressured 
lobe in the inner structure of the low-power radio galaxy Cen A. The closer 
such a structure is to a spherical expansion, the more normal the shock will 
be everywhere and the better the applicability of the above equations. 

5.2 Low-power FRI jets 

The appearance of low-power jets, in FRI radio galaxies, is quite different from 
that of high-power jets: compare Fig. with Fig. [HJ Whereas the high-power 
jets in FRII radio sources are generally weak features within the well-defined, 
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and much brighter radio lobes, low-power jets are usually of high contrast 
against the inner radio structures. Low-power jets are usually brightest near 
the active galactic nucleus, and then fade gradually in brightness at larger 
distances from the core, although this pattern is often interrupted by patterns 
of bright knots in all wavebands. FRI jets show a range of morphologies, from 
almost straight and symmetrical two-sided jets to the bent "head-tail" sources 
that are characteristic of low-power extended radio sources in clusters. 

While it is believed that the same basic mechanism is responsible for the 
generation of low-power and high-power jets, it appears that low-power jets 
dissipate much of their energy without developing a well-defined beam head. 
This is often taken to indicate that low-power jets are in good contact with the 
external medium in which they are embedded, so that they share momentum 
and energy with entrained material as they propagate. The strong velocity 
shear between a fast-moving (originally relativistic, based on VLBI observa- 
tions of the cores of low-power radio sources jSU]) jet flow and the almost 
stationary external medium must generate instabilities at the interface |10| . 
and drive the flow into a turbulent state. The physics of the resulting flow is far 
from clear, although with simplifying assumptions |Hj it has been shown that 
subsonic and transonic turbulent jet flows can reproduce the general trend 
of radio brightness and jet width for plausible confining atmospheres. More 
recently, model fits |S7] using (|41() and l|45l) to X-ray measurements of the 
thermal atmosphere of 3C 31, coupled with results from high-dynamic range 
radio mapping [ZH] (Fig. EJ, have led to a self-consistent model of the flow. 

The kinematic model of 3C 31 jZH] involves three distinct regions in the jet: 
an inner ~ 1 kpc long section of narrow jet where the flow is fast (relativistic) 
and the opening angle of the jet is small, a flaring region of ~ 2 kpc during 
which the jet broadens rapidly and brightens in the radio, and an outer region 
in which the jet expands steadily with a smaller opening angle than in the 
flaring region. Here the jet decelerates steadily from a moderate Mach number 
as it picks up mass from the external medium or from the mass lost by stars 
embedded in the jet [73j. Buoyancy forces are important for much of the 
outer flow, and the jet becomes sensitive to local changes in the density of the 
intracluster medium, and hence is liable to deflect from straight-line motion. 

While this model provides a good description of 3C 31, which displays a 
smoothly-varying radio brightness profile and so is amenable to simple fitting 
for flow speed, and is based simply on the basic conservation laws for mass, 
momentum, and energy |5| , it represents only the overall deceleration of the jet 
from entrainment. This leaves unanswered important questions about the ori- 
gins of the bright knots which appear in many jets, and which are interpreted 
as the sites of strong shocks (which are relativistic, based on proper motion 
studies) , and about the detailed physics responsible for mass entrainment and 
the deceleration of the jet. 

While the knots in jets are sites where intense acceleration of particles to 
highly relativistic energies can take place, it is clear from the diffuse X-ray and 
other emission between the knots that there is continuing particle acceleration 
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even in the inter-knot regions. The structures of the acceleration regions can 
only be studied in detail in the closest radio galaxies, such as Cen A jSS] 
and M 87 |!J5]- A number of features of the knots have become apparent 
through detailed multiwavelength studies, of which the most important are 
the presence of high-speed motions (apparent speeds greater than the speed of 
light having been noted in M 87) , and the distinct locations for the acceleration 
of particles radiating in different wavebands. 



6 The interpretation of multiwavelength data in terms of 
physical parameters 

6.1 Minimum energy and magnetic field strength 

The evidence for magnetic fields in the jets and lobes of active galaxies 
(Sect. Ifci.Bfl confirms the synchrotron origin of the radio emission, which is 
therefore an inseparable function of the magnetic-field strength and electron 
spectrum. To progress further it is usual to assume that the source is radiating 
such that its combined energy in particles and magnetic field is a minimum 
|20| . In this situation the energy in the magnetic field is ~ 3/4 of the energy 
in the particles, and so this is similar to the condition in which the two are 
equal and the source is in "equipartition" . 

We can calculate the minimum- energy magnetic field for a power-law spec- 
trum using the equations in Sect. 12. H and for more complicated spectra results 
can be obtained via numerical integration. The total energy in electrons 

/■7max Hm B x 

U c = Km c c 2 / 77~ p <i7 = k m c c 2 / j'^d-y (52) 

J 7 m i n « 7min 

From J7J), the synchrotron luminosity per unit frequency is 

L v (x nv- a B a+1 (53) 
and we can eliminate K from (|52|) and l|53|) to give (for a ^ 0.5) 

U e ex L^B-W (7^ - 7^„ 2 ") (1 - 2a)- 1 (54) 

If K is the ratio of energy in heavy non-radiating particles to that in 
electrons, and r\ is the filling factor in the emission region of volume V, the 
energy density in particles can be written as 

u p d = CiB ^ > ^ y L v v ( 7max -7 min M 1 - 2 ") (55) 

where the constant of proportionality, Ci, is determined from synchrotron 
theory and can be expressed explicitly using the equations in Sect. 12.11 The 
energy density in the magnetic field can be written as 
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u B = C 2 B 2 



(56) 



where the value of the constant C 2 depends on the system of units. The value 
of B which minimizes the sum of the two energy densities, B mc , is then given 

by 



LlyV 



-, l/(a+3) 



2C 2 r\V (l-2a) 
Since luminosity density and flux density S Vl are related by 



(57) 



L v = (1 + z) 0-1 S^nDl 



(58) 

where Dl is the luminosity density, and the volume of a radio source can be 
specified IHOJ in terms of its angular component sizes, £? x , 9 y and path length 
through the source, d, as 



V = 



dDl/{l + zf 



(59) 



we find 



O + iKMi + x) ^ 3+a (7,^ a "7^ n 2a ) n ' 

2C 2 " ( } (1 - 2a) 

(60) 

At minimum energy, the total energy density is u mc = u^{a + 3)/(o; + 1) = 
ttpci(o< + 3)/2. Any change in the ratio of ub to u pc \ increases the total energy 
and pressure of the source. Thus if, as seen in the southwest radio lobe of 
Cen A fSect. l6*3l . the minimum pressure is below that of the external medium, 
expansion can immediately be inferred if lobes are to stay inflated for an 
appreciable time. Using minimum energy, the magnetic fields in radio sources 
are generally estimated at 1-10 nT (10-100 /^Gauss). 

The dependence of B me on the lower and upper cut-off energies of the 
electron spectrum are weak, and since a is expected to be ~ 0.6 from shock 
acceleration p^, or larger where energy losses have steepened the spectrum, 
it is the lower energy which has most effect. There is a stronger dependence 
of B mc on (1 + K), the ratio of energy in all particles to that in radiating 
electrons, and r\. It is usual to adopt unity for both of these factors, which 
gives a true minimum value for the energy and pressure [(7 — 1) times the 
energy density, where in this context 7 is the ratio of specific heats, which 
is 4/3 for a relativistic gas]. Indirect arguments are then applied to estimate 
if (1 + K) should be larger than unity or rj less than unity for a particular 
source. 

In the case of SSC, L v must be converted into a spectral number of photons 
per unit volume, n v [i.e., n(v ) in the equations of Sect. |2~2"] . For a uniform 
spherical source, of radius r, the average value for n v is given by 
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4 c hv 

It is instructive to consider the effects of relativistic beaming (Sect. 01. For 
the case of a spherical blob, 1)58(1 is now the equation for L v 8^ +a \ and the 
total synchrotron luminosity in the source frame is reduced by a factor of S 4 , 
with ub and re, the normalization of the electron spectrum, both reduced by 
S 2 . SSC in the source frame is reduced by a factor of 5 6 , and thus sources which 
are more beamed have lower SSC to synchrotron flux ratios. [The factors of 
S 2 become 5 2 ( 2+Q )/( 3+a ) for the jet case where £( 2+a ) should be used.] 

The minimum-energy assumption can be tested by combining measure- 
ments of synchrotron and inverse-Compton emission from the same electron 
population. If the inverse Compton process is responsible for most of the 
higher-energy radiation that is measured, and the properties of the photon 
field are well known, the electron population is probed directly (Sect.J^J. In 
combination with measurements of radio synchrotron emission from the same 
electrons, both the electron density and magnetic field strength can then be 
estimated, and the minimum-energy assumption can be tested. Since the ap- 
plication of this test requires that the volume and any bulk motion of the 
emitting plasma be known, the best locations for testing equipartition are 
the radio hotspots, which are relatively bright and compact, and are thought 
to arise from sub-relativistic flows at jet termination (but see 0ED> and old 
radio lobes that should be in rough overall pressure balance with the external 
medium. There is no reason to expect dynamical structures to be at minimum 
energy. 

Chandra has allowed such tests to be made on a significant number of 
lobes and hotspots, with results generally finding magnetic field strengths 
within a factor of a few of their minimum energy (equipartition) values for 
(1 + if) = rj = 1 (e.g., [T2I ES|). A study of ~ 40 hotspot X-ray detections 
concludes that the most luminous hotspots tend to be in good agreement 
with minimum-energy magnetic fields, whereas in less-luminous sources the 
interpretation is complicated by an additional synchrotron component of X- 
ray emission |59| . 



6.2 Jet composition 

Much of the detailed physics of jets depends on what they contain. While it 
is clear from their polarized radiation that jets contain fast charged particles 
and magnetic fields, it is less clear whether these are primary energy-carrying 
constituents of the jet or secondary. A number of different possibilities for the 
energy carriers have been suggested. The most popular are that the jets are 
primarily composed of electrons and positrons, or of electrons and protons, 
although electromagnetic (Poynting flux) jets |1U8| and proton-dominated jets 
|MI I104| have been discussed. 
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Since jets are presumed to obtain much of their energy from the infall 
of matter into a supermassive black hole, it is natural to suppose that elec- 
tromagnetic radiation would carry much of the energy from the system on 
the smallest scales, since a plausible mechanism for the extraction of energy 
from the system is the twisting of magnetic field linked to the accretion disk. 
jH2| describe the processes that might launch such a jet, but fast interactions 
with the plasma environment and efficient particle acceleration are expected 
quickly to load the field with matter. In the resulting magnetohydrodynamic 
flow much of the momentum would be carried by particles, although Poynting 
flux may carry a significant fraction of the total energy | 5. . 

An electron-positron pair plasma is a natural consequence of the high 
energy density near the centres of active galaxies, and hence it might be 
expected that electrons and positrons would be an important, and perhaps 
dominant, component of the jet outflow. Acceleration of such a plasma by the 
strong radiation field near the active nucleus is certainly possible jS^l I1U9| , 
although radiation drag is an important limitation on the speeds to which 
the flow can be accelerated if only electron-positron pairs are present |121j . 
The importance of electron-positron pairs can be established by assuming 
that they are the major constituents of the flow, and then testing that the 
kinetic energy that they carry down the jet is comparable with the radiated 
power. While an electron-positron jet with arbitrary distribution of Lorentz 
factors, N("f), can satisfy this condition, an electron-proton jet should not 
extend to Lorentz factors 7 m in < 100, to avoid a high order of inefficiency 
and depolarization by Faraday rotation (Sect. Ifi.fi [I . Another test, based on 
the assumption that the hotspots at jet termination in the FRII radio galaxy 
Cygnus A are reverse shocks in the jet fluid [72], suggests that the jet material 
is an electron-positron fluid at that point. 

Further information supporting the interpretation of jets as electron- 
positron plasmas at their point of injection has been obtained by VLBI polar- 
ization studies |134| . Here the detection of circular polarization, interpreted 
as arising from mode conversion, strongly suggests that, on the pc scale, jets 
have low 7 m ; n (< 20) and hence are electron-positron plasmas (Sec. I6.6fl . al- 
though 110 find that other possibilities (notably that the jet is composed of 
electrons and protons) are not ruled out since other mechanisms exist for the 
generation of circular polarization. 

There are also grounds for expecting high-power jets to contain a consid- 
erable fraction of their mass in protons. Protons outnumber electrons by 100 
to 1 in Galactic cosmic rays at energies at which they leak out of the galaxy 
before significant energy loss |138| . Some jet acceleration mechanisms can ac- 
celerate heavy components of the jets as much as, or more than, the lighter 
components, and, even if the jet is initially light, interactions with the external 
medium are expected to load it with protons (Sec.J^J). However, there is little 
direct evidence for the proton contents of jets. Polarization studies do not find 
evidence of embedded thermal material, and protons are relatively inefficient 
radiators so the synchrotron emission of relativistic protons is probably not 
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detectable (but see 2 ). It does seem to be clear that jets are of low density 
relative to the external medium, on the basis of numerical simulations, since 
the structures of dense jets are unlike those observed, but there could still be 
appreciable mass in the flow. 

Support for the presence of protons is also found in models for the broad- 
band spectral energy distributions of emission from some sub-parsec scale 
quasar jets, where a significant proton contribution is required to boost the 
kinetic power sufficiently to match the total radiated power |12fi| . Some ra- 
dio lobes in low-power radio galaxies, if assumed to be radiating at minimum 
energy (Sect. [OJ, would collapse because the X-ray-emitting medium has a 
higher pressure Although there are several ways of boosting the internal 
pressure in such a situation, magnetic dominance would make the sources un- 
usual, electron dominance can be ruled out by constraints on inverse-Compton 
scattering of the CMB, and non-relativistic protons are disfavoured on grounds 
of Faraday rotation (Sect. leaving a relativistic proton component most 
likely. In contrast, there are other cases where the radio lobes would be over- 
pressured with respect to the ambient X-ray-emitting medium if a relativistic 
proton contribution were included jSE]. Although this in itself is not a major 
difficulty, since there are open issues concerning lobe expansion (Sect. I6.5JI . 
these are lobes for which the magnetic field can be measured (Sect. I6.1|) and 
the sources are close to minimum energy when only the radiating electrons 
are considered. 

It has been known for some time that the minimum pressure in low-power 
jets (calculated assuming an electron-positron plasma) is typically below that 
of the external X-ray-emitting medium, e.g., |421 1711 l9"Tl H4Uj . However, this 
cannot be used simply to infer that the jets are launched as an electron-proton 
plasma to give them the extra required pressure. Low-power jets are believed 
to slow down to sub-relativistic speeds via entrainment of thermal material 
(Sect.E2J, and even the most detailed hydrodynamical modeling, such as that 
which has been applied to 3C31 by does not decide the issue of primary 
jet content. 

On balance it appears that the dominant energy-carrying constituents of 
jets are electrons and positrons, but that relativistic protons may also be im- 
portant, particularly in momentum transport. However, the situation remains 
unclear, and there are degeneracies between the measurement of jet speed and 
jet composition which render this conclusion tentative. 

6.3 Jet speed 

The resolved X-ray jets (excluding hotspots) that Chandra detects in powerful 
radio sources are mostly in quasars, e.g., [2E1 IH3 [H2I EH EH EH EHI 
I12()j , with the bright radio sources Pictor A and Cygnus A [1361 I137| being 
exceptions. The quasar X-ray jet emission is one-sided, on the same side as 
the brighter radio jet, implying that relativistic beaming is important. Two- 
sided X-ray emission, such as that in the quasar 3C 9 most likely does 
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not imply the presence of counter-jet emission, but rather the presence of 
the more isotropic lobe, hotspot, or cluster-related emission expected at some 
level in all sources and detected in many, e.g., |56l I145J . 

Currently there are several tens of quasar X-ray jet detections. They have 
mostly been found through targeted programs to observe bright, prominent, 
one-sided radio jets. In most cases there has been no pre-existing reported 
optical jet detection, but there has been reasonable success from follow-up 
work. The level of many such detections lies below an interpolation between 
the radio and X-ray spectra, suggesting that synchrotron emission from a sin- 
gle power-law distribution of electrons is not responsible for all the emission, 
e.g., |113M117| . although it has been pointed out that since high-energy elec- 
trons lose energy less efficiently via inverse Compton scattering in the Klcin- 
Nishina regime, the electron spectrum should harden at high energies and in 
some cases may produce a synchrotron spectrum that can match observations 

m 

In order to avoid the total energy in particles and magnetic field being 
orders of magnitude above its minimum value, as would arise from a simple 
SSC explanation, the most widely favored model for the X-ray emission in 
these cases is that the X-rays are produced by inverse Compton scattering 
of CMB photons by the electrons in a fast jet (Sect. |3J that sees boosted 
CMB radiation and emits beamed X-rays in the observer's frame |27l I127| . 
The model can produce sufficient X-rays with the jet at minimum energy, but 
only if the bulk motion is highly relativistic (bulk Lorentz factor, r « 5 — 20) 
and the jet at a small angle to the line of sight. Although such a speed and 
angle are supported on the small scale by VLBI measurements, at least for 
the source which has guided this work, PKS 0637-752 H3HH3, the jet must 
remain highly relativistic hundreds of kpc from the core (after projection 
is taken into account) for the X-rays to be produced by this mechanism. 
This conclusion, based on multiwavelength data, has been something of a 
surprise, since earlier statistical studies of the structures of powerful radio 
sources suggested that jet velocities average only about 0.7c at distances of 
tens of kpc from the core I133j . 

A difficulty with the interpretation of the X-rays from radio jets as due to 
inverse-Compton scattering of CMB photons by a fast jet is that the observed 
gradients in X-ray surface brightness at the edge of the knots are sharp |128j . 
Since the X-rays are generated from low-energy electrons (as the CMB photons 
are boosted in the electron rest frame), the lifetimes of the electrons are long 
and it is difficult to see how X-ray knots in high-power jets can have a steeper 
gradient (as in PKS 0637-752, |117j ) than radio knots, which are produced by 
electrons of similar or higher energy. Although |128j suggest that the effects of 
strong clumping in the jets may resolve this issue, a fast jet and the proposed 
mechanism is then no longer required, since such clumping would increase the 
SSC yield for a slow jet at minimum energy |117| . 
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6.4 Particle acceleration in jets 

The electrons that emit synchrotron radiation in the radio, optical, and X-ray 
wavebands are all ultrarelativistic (with Lorentz factors > 10 3 ). The existence 
of large numbers of such relativistic electrons depends on their acceleration to 
high energies locally within the jet, or hot-spot, since their lifetimes against 
synchrotron losses (Sec. I2.1f> are usually less than the minimum transport 
times from the active nuclei. (This may not be the case if proton synchrotron 
radiation is important 0.) Particle acceleration is generally discussed for the 
cases of a particle interacting with a distributed population of plasma waves 
or magnetohydrodynamic turbulence, or shock acceleration. Reviews of these 
processes may be found in ^3 [SHI E3 and elsewhere. 

Resolved X-ray jets in active galaxies with low radio power are detected 
with Chandra in sources covering the whole range of orientation suggested 
by unified schemes, suggesting that beaming is less important than in their 
more powerful counterparts. The several tens of detected sources range from 
beamed jets in BL Lac objects 1101] to two-sided jets in radio galaxies 
P91 135] . with most X-ray jets corresponding to the brighter radio jet, e.g., 
[SSI E3 EOl EH ESI EHI Several of the observations have been targeted at 
sources already known to have optical jets, from ground-based work or HST. 
However, it's easier |144| to detect X-ray jets in modest Chandra exposures 
than to detect optical jets in HST snapshot surveys, because there is generally 
better contrast with galaxy emission in the X-ray band than in the optical. 

Inverse Compton models for any reasonable photon field suggest an un- 
comfortably large departure from a minimum-energy magnetic field in most 
low-power X-ray jets, e.g., |55) . Synchrotron emission from a single electron 
population, usually with a broken power law, is the model of choice to fit the 
radio, optical, and X-ray flux densities and the relatively steep X-ray spectra, 
e.g., |16U55| . X-ray synchrotron emission requires TeV-energy electrons which 
lose energy so fast that they must be accelerated in situ. 

The above arguments, applied to the bright northeast jet of the nearest 
radio galaxy, Cen A, find in favor of X-ray synchrotron emission |74| . and the 
proximity of Cen A allows its acceleration sites to be probed in the great- 
est possible detail. Unfortunately the dramatic dust lane spanning the galaxy 
masks any optical jet emission. Proper motion of order 0.5c, observed both 
in the diffuse emission and some knots of Cen A's radio jet, is indicative of 
bulk motion rather than pattern speed [SB]- Since Cen A has a strong jet 
to counter-jet asymmetry, Ij23(l then suggests that the jet is at a small angle 
to the line of sight. Since this contradicts the evidence based on parsec-scale 
properties, and other considerations, that the jet is at about 50 degrees to the 
line of sight, Cen A appears to be a case where intrinsic effects render 123|) 
inapplicable. Some of the bright X-ray knots have only weak radio emission 
with no indication of proper motion, but with the radio emission brightening 
down the jet in the direction away from the nucleus. While the radio asso- 
ciation confirms that these X-ray knots are indeed jet related, the emission 



26 D.M. Worrall & M. Birkinshaw 



profiles are not what are expected from a simple toy model where the elec- 
trons are accelerated and then advect down the jet, with the X-ray emitting 
electrons losing energy faster than the radio-emitting electrons. Instead, it is 
proposed [HE] that there are obstacles in the jet (gas clouds or high-mass-loss 
stars). Both radio and X-ray-emitting electrons are accelerated in the stand- 
ing shock of this obstacle, and a wake downstream causes further acceleration 
of the low-energy, radio-emitting, electrons. The resulting radio- X-ray offsets, 
averaged over several knots, could give the radio- X-ray offsets seen in more 
distant jets, e.g., [53] , 

The X-ray and radio emission in hot-spots are also offset in some cases 
|5B"] , presumably also because of the different locations of acceleration of the 
particles at the different energies probed in these wavebands, if the emission 
is all synchrotron. However, it is expected that much of the X-ray emission 
should be of inverse-Compton origin, and then the offsets become harder to 
understand, although have described a model in which decelerations of 
the jet plasma near a hot-spot can generate X-ray enhancements and small 
offsets. 

Optical polarization might help us to learn more about the acceleration 
processes in low-power jets. In M 87 there is evidence for strong shock ac- 
celeration at the base of bright emitting regions, in compressed transverse 
magnetic fields A knot in the jet of M 87 has been observed to vary 
in the X-ray and optical on the timescale of months, consistent with shock 
acceleration, expansion, and energy losses |62l 1100) . 

6.5 Pressure and confinement 

It is interesting to compare the minimum pressure (Sect. 16. 1[) in radio lobes 
with that of the external X-ray-emitting environment. Over-pressure in the 
lobe would imply an expansion, which may be supersonic and should involve 
significant heating of the external gas. Under-pressure suggests either that 
the lobe is undergoing collapse (this should be rare given their prevalence) or 
that there is an additional component of pressure that may be in the form of 
relativistic protons. FRII sources are generally at high redshift where obser- 
vations lack sufficient sensitivity and resolution to draw strong conclusions. 
However, notwithstanding the fact that there is no reason to expect dynam- 
ical structures to be at minimum energy, where tests are possible it tends 
to be confirmed fSect. l6~T)l . and it appears that with this assumption rough 
pressure balance prevails, e.g., [71 1331 . 

The situation in low-power radio galaxies is more complicated, because the 
medium plays an important role in the deceleration of the jets, such that they 
share momentum and energy with entrained material. However, these sources 
have the advantage of being closer, and Sect. 15.21 describes how the pressure 
profile predicted by the mass-entrainmcnt model for deceleration of the jet in 
3C 31 gives an excellent match to observations. The outer radio structures of 
FRIs may sometimes be buoyant, e.g., |143) . and in other cases show evidence 
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of having done significant work on the gas in pushing it aside |151 I89| , or 
responsibility for lifting gas in hot bubbles, e.g., [51]. expected to result in 
eventual heating, e.g., |ll)5| . Such heating would help to explain the weakness 
or absence of lines from gas cooling below 1 keV in the densest central regions 
of galaxy and cluster atmospheres, e.g., jl02j . A statistical study shows that 
atmospheres containing radio sources tend to be hotter than those without 

ma 

One place where heating is definitely expected is from gas crossing the 
supersonically advancing bow shock of an expanding lobe (Sec. 15.1(1 . It is 
possible to interpret X-ray cavities coincident with the inner parts of the radio 
lobes of Cygnus A as due to an emissivity contrast between bow-shock heated 
gas outside the lobes and the more easily detected ambient cluster medium 
|24| . although the parameters of the shock are not effectively constrained by 
the data. More recent Chandra observations of Cygnus A find gas at the sides 
of the lobes to have kT ~ 6 keV, slightly hotter than the value of 5 keV 
from ambient medium at the same cluster radius, possibly indicating cooling 
after bow-shock heating, but again the data do not usefully constrain model 
parameters |122j . 

The first and best example of a shell of heated gas which can reasonably 
be attributed to supersonic expansion is not in an FRII radio source, but in 
Cen A (Fig. El). High-quality Chandra and XMM-Newton data [7J] provide the 
temperature and density constraints needed to test the model and measure 
the supersonic advance speed of the bow shock responsible for the heating. 
This source is an excellent example where much of the theory outlined earlier 
in this chapter can be applied. 

Cen A is our nearest radio galaxy, at a distance of 3.4 Mpc [S7J so that 
1 arcsec corresponds to ~ 17 pc. The full extent of Cen A's radio emission 
covers several degrees on the sky (HH]- Within this lies a sub-galaxy-sized 
double-lobed inner structure [2] with a predominantly one-sided jet to the 
northeast and a weak counter-jet to the southwest [SS], embedded in a radio 
lobe with pressure 1.4 x 10~ 12 Pa or more, greater than the pressure in the 
ambient ISM (~ 1.8 x 10~ 13 Pa; Tabled]), and so which should be surrounded 
by a shock. Around this southwest lobe there is a shell of X-ray emitting gas 
which appears to have the geometry of the shocked ambient gas in Fig. H3 
[25]. Although the capped lobe is around the weak counterjet, so it is not 
evident that the lobe is being thrust forward supersonically with respect to 
the external interstellar medium (ISM) by the momentum flux of an active 
jet, the current high internal pressure in the radio lobe ensures its strong 
expansion. 

The temperature, proton density and pressure of the ambient ISM and 
the X-ray shell, taken from [73] are given in Table The ambient medium 
is measured to have n p ~ 1.7 x 10 3 m -3 and kT — 0.29 keV, whereas the 
shell is ten times hotter, at kT — 2.9 keV, and twelve times denser, with 
n p ~ 2 x 10 4 m -3 . From f5T)l and lf50"|l . we see that temperature and density 
measurements for both the ambient medium and the shocked gas directly test 
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Fig. 9. 0.5-5 keV X-ray image of Cen A from a 50 ks Chandra exposure using 
the ACIS-S instrument. Pixel size is 1 arcsec. Data are exposure-corrected. Many 
point sources associated with Cen A's host galaxy NGC 5128 are seen, together 
with a diffuse background of X-ray-emitting gas of kT » 0.3 keV, and labelled X- 
ray structures that are in whole or part related to radio structures (see the text and 
EHlESHHIZg for more details) 

shock heating, since only two of the four parameters are required to measure 
the Mach number, and the other two test the model. 

The most straightforward application of the equations finds that the densi- 
ties and temperatures are not self-consistent. The shell's density and tempera- 
ture are wrong for gas directly in contact with the bow shock. However, we can 
find a Mach number consistent with shocking the gas to a temperature and 
density such that the combined thermal and ram pressure is in pressure equi- 
librium with the thermal pressure of the detected shell: M. = 8.5, w a dv ~ 2400 
km s _1 . The post-shock temperature is kT^ ~ 6.8 keV. The 6.8 keV gas flows 
back from the shock, into the X-ray-detected shell at 2.9 keV. The character- 
istics of this undetected hotter gas are given in Tabled I n this table we also 
quote estimates of the relative X-ray emissivity (per unit volume) of gas in 
the different structures over the 0.4-2 keV energy band, where the Chandra 



Multiwavelength evidence of the physical processes in radio jets 29 
Table 1. Physical parameters of the gas in various regions of Cen A 



Structure kT n p Pressure 0.4-2 keV 

(keV) (m -3 ) (Pa)f relative 

emissivity, £ 



ISM (measured) 


0.29 


1700 


1.8 x 10 _13 i 


1 


Behind bow shock (inferred) 


6.8 


6530 


2.1 x 10~ n 


13 


Shell (measured) 


2.9 


20000 


2.1 x 10" 11 


127 



f 1 Pascal = 10 dyn cm 2 J incorrectly reported in Table 5 of 75 

ISM and Shell pressures are thermal only. Inferred pressure includes ram pressure 
ofpiv 2 dv (M 2 + 3)/4.M 2 

The minimum-energy pressure in the radio lobe is ~ 1.4 x 10 12 Pa. 



response peaks and is relatively flat. The gas directly behind the bow shock 
has a predicted emissivity that is an order of magnitude fainter than that in 
the shell, accounting for its non-detection in our measurements. 

The radiative timescale for material in the shell is ~ 2 x 10 9 yrs (Fig. [SJ, 
which is large compared with the lobe expansion time (< 2.4 x 10 6 years), 
so the material in the shell behaves as an adiabatic gas 3 . The shell is over- 
pressured compared with the minimum-energy pressure in the radio lobe (in 
magnetic field and radiating electrons) by a factor of ~ 10. If we assume min- 
imum energy in the lobe, despite it being a dynamical structure, and that the 
shell has reached equilibrium [but note that the sound-crossing time in the 
shell (thickness ~ 0.3 kpc, c s ~9x 10~ 7 kpc yr _1 ) is about 15 per cent of the 
maximum time we estimate it has taken the lobe to reach its current size] , the 
shell's overpressure relative to the radio lobe could be balanced by the ram 
pressure from internal motions in the lobe for a moderate relativistic proton 
loading. Finally, the shell's kinetic energy is ~ 5 times its thermal energy, 
and exceeds the thermal energy of the ISM within 15 kpc of the centre of the 
galaxy. As the shell dissipates, most of the kinetic energy should ultimately 
be converted into heat and this will have a major effect on Cen A's ISM, 
providing distributed heating. 



6.6 Magnetic field structures 

Jets contain significant magnetic fields. The orientation of these fields is dis- 
played by the linear polarization that they show in their radio synchrotron 
emission, and for many years detailed maps of the radio polarizations of jets 
have been used to infer the magnetic field geometries both on kpc scales, e.g., 
for NGC 6251 [HE| and 3C 31 g3|, and pc scales, e.g., for BL Lac objects 

The fractional linear polarization expected from an optically-thin synchro- 
tron-emitting plasma is 

a + 1 . . 

= — T 62 
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where a is the radio spectral index, and it is assumed that the frequency 
of observation satisfies J^x n v e <C v <C Tmax^g) as f° r 0- The electric field 
of the emitted radiation is perpendicular to the magnetic field within which 
the emitting electron population lies, so that an observation of the plane of 
polarization of the radio emission can be used to infer the geometry of the 
magnetic field within the radio jet. 

In mapping the direction of the projected magnetic field, account has to be 
taken of the effect of Faraday rotation of the plane of polarization by mixed 
thermal material and fields within and around the radio jets. For a foreground 
plane slab of plasma, at a wavelength of A, the Faraday rotation is an angle 

F = RM A 2 (63) 
where the rotation measure along path length L is 

1 f L 

RM = - / n 0jth r c u g cos-0 dz (64) 
c Jo 

where n c th is the thermal electron density, r c is the classical electron radius, 
z/ g is the non-relativistic gyrofrcqucncy, and ?/> is the angle of the magnetic 
field to the line of sight. In SI units, 

RM = 8100 [ u e ,th B\\ dz rad m~ 2 (65) 
Jo 

for a magnetic field with line-of-sight component Bii (in Tesla), where n C! th 
is in m~ 3 and the thickness L is in pc. If the thermal material is mixed with 
the radio jet's plasma, an effect half this size is expected. If the magnetic field 
and thermal material have a complicated structure, different dependencies of 
9-p on A are possible. Faraday rotation is measured by multi-frequency radio 
mapping, and uncertainties in the intrinsic plane of polarization of the radio 
emission are minimized by working at high radio frequencies. 

Since the presence of an appreciable density of thermal material around the 
jets is revealed in X-rays (Sect. UJ, a significant Faraday rotation is possible 
if this medium contains a magnetic field, and some evidence for such Fara- 
day rotation has been seen, e.g., in Hercules A [63], with the brighter, and 
presumably Lorentz-boosted, jet often showing a lower Faraday rotation than 
the opposite jet, as might be expected from simple geometrical considerations 

ma E7]. 

A more important consequence of the significant linear polarizations mea- 
sured for radio jets, however, is that those jets cannot contain large amounts 
of thermal material. If they did, then the change in Of between the front and 
back sides of the jets would lead to linear polarizations far less than the ob- 
served values, which are sometimes near the theoretical maximum given by 
<|62|) . The limits on Faraday depolarization thus lead us to conclude that the 
X-rays from radio jets and lobes cannot come from thermal material mixed 
with the jets. 
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Polarization mapping also reveals that the magnetic fields in radio jets are 
relatively well ordered. In many jets of both high and low radio power the 
magnetic fields are longitudinal for the first few kpc of the length of the jet, 
then, in lower-power jets, the magnetic field becomes perpendicular to the 
jet axis at larger distances from the core, often after a significant change in 
the jet's width [17) . Variations on this pattern are seen where bright emission 
knots exist in the jet. These knots appear to be associated with strong shocks, 
which compress the magnetic field strongly into a transverse pattern at the 
upstream (core-side) of the knot, after which the magnetic field pattern be- 
comes complicated. A particularly good example is M 87, where the use 
of HST to map the optical polarization provides additional information, since 
the lifetimes of the electrons emitting optical synchrotron emission are far 
shorter than those of radio-emitting electrons. Significant differences between 
the polarization structures seen in the optical and radio in M 87 suggest that 
the sites of acceleration are different for different electron energies, with the 
strongest shocks, that provide acceleration to the highest energies, appearing 
in the most central parts of the jet. 

In the central parts of radio sources, VLBI observations have been able 
to detect both linear and circular polarization. A distinct difference is seen 
in the polarizations of the pc-scale (VLBI) jets of BL Lac objects and radio 
galaxies or quasars. Generally, the magnetic field orientation in the cores of the 
highest-power sources is along the jet but in lower-power sources, such as 
BL Lac objects, regions of significant transverse field are seen, and change on 
short timescales |UJ- Again this is interpreted as the effect of shock structures 
moving along the jets, compressing the magnetic field as they pass, although 
helical patterns in the flow may affect the apparent field pattern, as may the 
effects of relativistic aberration. The lower linear polarizations often seen in 
pc-scale jets may be due to depolarization caused by the superposition of many 
small-scale structures, with different field orientations, within the resolution 
of the observations, or perhaps to significant structure in the nuclear gas near 
the jet. 

Circular polarization is generally undetectable in large-scale jets, but has 
been mapped in pc-scale jets using VLBI [""""] [""""J, where it is interpreted 
as arising from mode conversion — the conversion of linear polarization to 
circular polarization by Faraday rotation in the jet plasma, although this is 
not certain [110) . If the X-ray emission from powerful jets on the large scale is 
interpreted as inverse-Compton scattering of CMB photons in a fast jet, then 
the population of relativistic electrons has to extend down to low energies 
(Lorentz factors 7 m i n = 10 — 20, [2*71 1127) ). On the other hand, the detection 
of significant polarization in VLBI jets implies that there should be few low- 
energy protons and electrons (which would cause excessive internal Faraday 
rotation), and hence that the jets should be composed of an electron-positron 
plasma. However, the presence of an electron-positron plasma in the VLBI 
jet does not necessarily imply its presence on the largest scales, where the jet 
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may have become loaded with material gathered from the ambient medium 
or the radio lobe. 

6.7 Cores: the inner jets 

Radio jets extend into the cores of active galaxies, where they are faster and 
more compact. Special-relativistic effects then cause their brightness and vari- 
ability time scales to be strong functions of jet orientation. As a result of 
synchrotron self-absorption, it is important to use VLBI techniques at high 
radio or mm frequencies to see close to the bases of jets. 

The inner jets are difficult to distinguish at non-radio wavelengths because 
of their closeness to the central engine, which is bright at infrared to X-ray 
energies, and because of orientation-dependent absorption in the optical to 
soft-X-ray bands from gas and dust structures. These gas and dust structures 
are only sometimes detected, but their presence is often inferred so that models 
which unify various classes of active galaxy might work, e.g., They affect 
the central engine the most. Thus, the first problem in using multiwavelength 
data to gain additional insight into radio jets on small scales is to separate out 
the jet emission. For no waveband does telescope spatial resolution match that 
of radio VLBI, and so only the tool of spectral separation is usually available. 

In the most core-dominated radio sources, such as bright variable BL Lac 
objects and quasars (often classed together as blazars), jet emission appears 
dominant at all energies, sometimes up to TeV. The multi- wavelength spec- 
tral energy distributions and variability time scales are used to probe the 
beaming parameters and the physical properties of the emitting regions, e.g., 
[481 I7fil 1125) . Correlated flares are sometimes measured across wavebands, 
giving support to the present of a dominant spatial region of emission, e.g., 
[1241 1131] . but otherwise uncertainties of size scales, geometries, and parame- 
ters for the competing processes of energy loss and acceleration, often force the 
adoption of oversimplified or poorly-constrained models for individual jets. 

In the quasar population in general, there is good evidence that in the 
X-ray an inner jet dominates the emission from core-dominated quasars, but 
not lobe-dominated quasars, e.g., [181 171)1 11411 11421 I148| . Radio galaxies are of 
particular interest, because here emission from the central engine is weak. In 
these sources a correlation of core soft- X-ray and radio emission, lost in lobe- 
dominated quasars, re-appears |22EniE2IISni, suggesting that jet-related X- 
ray emission is again dominant. In low-power radio galaxies, an optical core is 
often seen with HST, and is interpreted as synchrotron emission from a similar 
small-scale emitting region 23, 30. 53l ll32| . However, at higher X-ray energies, 
a number of radio galaxies show clear evidence of a hard, highly absorbed 
continuum, sometimes accompanied by Fe-line emission, e.g., [511 H2& . Here 
the decreasing jet output is leading to the central engine becoming increasingly 
dominant. Both X-ray components can sometimes be distinguished in the 
same spectrum, e.g., |3lfll3"5] . 
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The richness of high-energy structure in larger-scale radio jets has been 
revealed because it is well resolved with HST and Chandra. Our knowledge 
of the inner jets is limited by the confusion of components that has now been 
lifted for the jets discussed earlier in this chapter. 

7 Conclusion 

The new results on radio jets which have resulted from complementary X-ray 
and optical observations have brought some surprises. Firstly, synchrotron X- 
ray jets are common in low-power sources, which implies that the intrinsic 
electron spectrum continues to TeV energies, and requires substantial in situ 
particle acceleration. Secondly, the detection of many quasar X-ray jets, most 
commonly interpreted as due to beamed CMB photons, would suggest that 
highly relativistic bulk flows exist far from the cores, contradicting earlier 
statistical studies of radio sources (541 I133| . Jet theory has had some pleas- 
ing successes, such as the agreement of the X-ray pressure profile with the 
prediction from a hydrodynamical model for 3C 31 |79| . 

There is still much observational work to be done. Firstly, there is consid- 
erable bias in the jets which have been observed in the X-ray, and we need 
observations of unbiased samples over broader luminosity and redshift ranges, 
together with observations of the X-ray-emitting medium through which the 
jets propagate. The measurements should elucidate the relative importance of 
synchrotron emission, inverse Compton scattering, and relativistic beaming, 
refine our knowledge of the source energetics, and improve constraints on jet 
composition and speed. Secondly, we need more deep X-ray observations (and 
refined theory) to understand jet-lobe/intercluster medium interactions. Fi- 
nally, to study acceleration sites and processes, deeper and more detailed mul- 
tiwavelength mapping, spectroscopy, and temporal monitoring is required. In 
combination with multifrequency polarization measurements, such data could 
map the spatial distributions and follow the acceleration of the electrons re- 
sponsible for the radiation in the radio to X-ray bands. 

It should, however, be acknowledged that much basic physics of jets is 
still not well understood, e.g., the origin of the magnetic field, the method of 
jet production and collimation, and the effects of turbulence. Much theoreti- 
cal work in these areas is necessary to relate measurements to the processes 
occurring in the fascinating jet flows that are observed. 
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